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Abstract - This paper presents an approach for reducing speckle 
in ultrasound images using Coupled Partial Differential 
Equation (CPDE) which has been obtained by uniting second- 
order and the fourth- order partial differential equations. Using 
PDE to reduce the speckle is the noise- smoothing methods 
which is getting attention widely, because PDE can keep the 
edge well when it reduces the noise. We also introduced a 
median regulator to guide energy source to boost the features 
in the image and regularize the diffusion. The proposed 
method is tested in both simulated and real medical 
ultrasound images. The proposed method is compared with 
SRAD, Perona Malik diffusion and Non linear coherent 
diffusion methods, our method gives better result in terms of 
CNR, SSEVI and FOM. 

Index Terms - Ultrasound, Speckle, PDE, median filter. 

I. INTRODUCTION 

ULTRASOUND speckle is the result of the diffuse 
scattering, which occurs when an ultrasound pulse randomly 
interferes with the small particles or objects on a scale 
comparable to the sound wavelength. Speckle is an inherent 
property of an ultrasound image, and is modeled as spatial 
correlated multiplicative noise. In most cases, it is considered 
a contaminating factor that severely degrades image quality. 
To improve clinical diagnosis, speckle reduction is generally 
used for two applications: visualization enhancement and 
auto-segmentation improvement. Most speckle filters are 
developed for enhancing visualization of speckle images [1]- 
[3] including temporal averaging, median filtering, adaptive 
speckle reduction (ASR) and wavelet based methods [4] -[5], 
frequency compounding methods [6]-[7]. However, above 
methods could not succeed to balance between speckle 
suppression and feature preservation due to the complexity 
of speckle statistics. The comparative analysis [8] and review 
methods [9] also says that a technique that relies on a more 
accurate model for removing speckle noise while preserving 
image feature well would be rather valuable for practical use 
improve image simplification, which is in turn very beneficial 
in automated object detection. 

II. BACKGROUND OF PDE 

The use of partial differential equations (PDE) in image 
processing has grown significantly over the past years. Its 
basic idea is to deform an image, a curve or a surface in a 
partial differential equation framework, and to approach the 
expected result as a solution to this equation. P. Perona and 
J. Malik put forward an anisotropic diffusion (AD) equation 
to smooth a noisy image: 
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du(x, yj) 

dt 



= div(g\Vu(x, y,t)\)Vu(x, y,t) 



(1) 



where mfcyX) :ax-<0 ,+«)-►* i s a sca i e image, g(IVul) 
is a decreasing function of the gradient. Their work made an 
important influence on this field. 2D SRAD takes the format 
of the PDE of conventional anisotropic diffusion. Given an 
intensity image I Q (x, y) having none zero- valued intensities 
over the image domain U, the continuous form of SRAD is 
expressed as follows [10] 



dl(x,y,t) 



dt 



= div [c(x, y, : t)VI(x, y, t)] 



(2) 



I(x,y:0) = I (x,yUdI(x,y;t)/8Q = 

where 5Q denotes the border of Q._ n^ is the outer normal to 
the dQ, c(x) is the diffusion coefficient, and q is the ICOV. 
Both P-M diffusion and SRAD introduce blocky effects in 
image. SoYu-Li andM.Kaveh [11] proposed a fourth-order 
partial differential equation. This equation can reduce the 
noise, at the same time it can keep the image edge better. The 
PDE is: 



di_ 

dt 



= -V 2 [g(V 2 /)V 2 /] = -VU(A/)] (3) 



Where 



g(x) = 



1 



H + (s/k) 2 ] 



, k is constant. 



The advantage of using fourth-order PDE in denoising image 
is it removes the blocky effects that made by second-order 
nonlinearity diffusion equation. 

EI. PROPOSED METHOD 

The coupled PDE based on second-order and fourth-order 
PDE is: 



di_ 

dt 



= -V[ gl (M)VI + a.g 2 (M)AVI] (4) 



Where 'a' is the coupled coefficient and the functions g 1? g 2 
are plus and non-increasing function. The median filter is a 
well-known "edge preserving" nonlinear filter. It removes the 
extreme data while producing a smoothed output. We propose 
to use a median regulator in coupled PDE, which combines 
both medina filtering and strength of unity of second order 
and fourth order Partial Differential equation. The proposed 
technique is defined by the following relations: 
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dl(x,y,t) 
8t 



= -V[ gl (AI)VI + 



a.g 2 (AI)AVI] + af 



i(x, y,t) t=0 = I 

d n I = , / = median (I). 



(5) 



Speckle noise is signal-dependent noise. Typically, the bright 
regions have stronger noise than the dark regions. The source 
term/provides a boosting force to guide (or normalize) the 
diffusion evolution and also accelerate the convergence rate 
compared to natural diffusion, a is control coefficient whose 
value lies between and 1 
The function g v g 2 are defined as 



g l (x) = exp(-[x/k l ] 2 ), 



8 2 (x) = 



1 



(l + (x/k 2 ) 2 ) 



(6) 



Where k , k 2 are constants. The coupled PDE removes the 
noise and keeps the edges, regulator act as an accelerating 
force to produce an image with less noise and enhanced 
structure 

IV. Implementation and Result 

The proposed method is tested on a simulated fetus image 
using Field-II ultrasound simulation software and a real 
ultrasound liver image. Throughout this paper, we take 
multiplicative noise as the model for speckle. The degraded 
intensity for each pixel will be 



h,=f, 



u ■ 

i, j i, j 



(7) 



Where f t j is the intensity at pixel ( i, j )in the original 

image, u tj the multiplicative noise, and I t j is the intensity 

at pixel () in the observed image. The structure similarity 
(SSIM), Contrast to noise ratio (CNR), Pratt's Figure of Merit 
(FOM) are calculated to compare the performance of proposed 
method with SRAD, Perona Malik diffusion and Non linear 
coherent diffusion methods. The structure similarity (SSIM) 
[12] was used to evaluate the overall processing quality 



SSIM= P/W+fiX^cJ 



cs) 



Where jU t and (J. are the mean and standard deviation of 

intensities of pixels in a local window, and is the constant to 
avoid instability. All involved parameters were set as 
suggested in [13]. The value of SSIM varies between andl, 
with unity for ideal processing quality, i.e., when the filtered 
image is equal to the reference image. The contrast-to-noise- 
ratio (CNR) or lesion signal-to-noise ratio [13], [14], was 
computed by 
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CNR = 



L"i-ft 



>R 



P) 



-hCT-n 



Where /i x and <j\ are the mean and variance of intensities 

ofpixels in a region of interest (ROI), and ju 2 and <j\ are the 

mean and variance of intensities of pixels in a background 
region that has the same size as the ROI to be compared with. 
Edge preservation or edge enhancement can be quantitatively 
evaluated by Pratt's figure of merit (FOM). 



FOM = 



max 






(10) 



Where N,N ideal are the numbers of detected and ideal edge 

pixels, respectively, di is the Euclidean distance between the 
ith detected edge pixel and the nearest ideal edge pixel. A is 
a constant typically set to 1/9. The dynamic range of FOM is 
between [0, 1]. Higher value indicates better edge matching 
between processed image and the ideal image. In the 
implementation of median regularized coupled PDE method 
the following values are selected for the parameters: ^=0.03, 
k 2 =0.04 and a=l .The mask size of 3X3 is chosen to calculate 
median regulator. The result of proposed method, NCD, 
SRAD and Perona Malik diffusion on Field II ultrasound 
simulation software generated fetus image are shown in fig. 
1 . The same methods are tested on and real ultrasound liver 
images and the results are shown in fig.2.Comparision of 
performance measure are shown in Table. 1. Image profile 
along 1 54 th row of various filtered images are shown in the 
fig.3 

V. DISCUSSION AND CONCLUSION 

The result processed by the new method is much sharper 
in terms of edge preservation and smoother in terms of 
speckle noise reduction than the other three filtered results. 
The execution time is also much shorter than the other 
methods. For quantitative quality evaluation, we provide 
three metrics SSIM, CNR and FOM. Table. 1. shows the 
evaluation results for the processed image. The FOM value 
indicates that the new method is better than other methods 
in terms of edge preserving ability. CNR and SSIM values 
indicate that the new method gives a better processing result 
in terms of structure preservation and contrast enhancement. 
In order to keep the image feature better we adopted coupled 
PDE by uniting second order and fourth order PDE, Further 
more we added a median regulator to accelerate the 
convergence rate. The proposed Median regularized coupled 
PDE provides improved result over other methods in terms 
many parameters. Therefore it has a wide applicability in 
practice. 
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Fig.l Result on simulated image (a) Original Noisy Image , Image filtered by (b) Perona Malik Diffusion 
(c) Non Linear Coherent Diffusion (NCD) (d) SRAD (e) Medial Regularized Coupled PDE 

Table 1 Comparison of Performance Measure 
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Fig. 2. Result on real ultrasound image (a) Original Noisy Image , Image filtered by (b) Perona Malik Diffusion 
(c) Non Linear Coherent Diffusion (NCD) (d) SRAD (e) Medial Regularized Coupled PDE 
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Fig. 3 Profile of the filtered image along 154 th row (a) Perona Malik Diffusion (b) Non Linear Coherent Diffusion (NCD) (c) SRAD 

(d) Proposed Medial Regularized Coupled PDE 
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